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摘 要 : 针对 树木 点 云 拓扑 
为 了 在 点 
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云 表面 直接 应 用 网 格 收缩 算法 ， 
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进行 局 部 主 成 分 分 析 和 Delaunay 三 角 剖 分 ;其 次 ， 针 对 树木 点 云 拓 


题 ， 用 则 率 法 线 流 算 子 对 点 去 进行 收缩 针对 树木 枝条 细 长 有 曲 幅 度 平缓 等 特点 ， 


利用 \ 改 进 后 的 QEM 网 格 简化 方法 将 三 角 网 格 折 县 成 一 维 曲线 骨架 ;最 后 ， 将 得 到 的 曲线 骨架 进行 连通 和 居中 处 理 。 


提出 的 算法 直接 在 点 云 上 进 


行 操作 ， 不 需要 额外 的 信息 和 预 处 理 操作 ， 
明 ， 该 算法 提取 的 树木 点 云 骨架 充分 表达 了 树木 在 自然 环境 下 的 生物 性 结构 和 特征 ， 相 对 于 rosa、L1- 中 轴 等 经 
法 ， 在 树木 点 云 的 骨架 oe ne 枝条 重 
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时 骂 声 和 残缺 点 云 有 良好 的 鲁 棒 性 。 实 验证 


典 算 


Abstract: Aiming at that complex topological structure and various feature details of tree point cloud, this paper proposes 


an algorithm for extract the curve skeleton based on point cloud contraction. First, In order to directly apply the mesh 


shrinkage algorithm on the surface of the point cloud, performing local cloud principal component analysis and Delaunay 


triangulation on the point cloud. Secondly, for the problem that the tree point cloud’s complex topology and the details of 


the last branch, the curvature normal operator is used to shrink the point cloud. In view of the slenderness of the branches of 


the trees and the gentle curvature, the modified QEM mesh simplification method is used to fold the triangular mesh into a 


one-dimensional curve skeleton. Finally, connecting and centering the resulting curve skeleton. The algorithm in this paper 


operates directly on the point cloud and does not require additional information and pre-processing operations. Good 


robustness to noise and residual fault clouds. Experiments Show that compared with other classical algorithms such as L1 


and rosa, the tree point cloud skeleton extracted by the algorithm has a good topological structure, which fully expresses the 


biological structure and characteristics of trees in the natural environment. The skeleton extraction speed of the tree point 


cloud is increased by more than 3 times, and the branch reconstruction degree is increased by 25%. 
Key words: curvature normal flow operator; point cloud contraction; PCA; Delaunay triangulation; QEM 
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据 进 行 全 局 收缩 。 在 每 次 迁 代 收缩 中 ， 用 不 同 权重 


环 邻 域 边界 权重 ， 
在 处 理 不 完整 点 云 和 了 品 
采样 半径 的 大 小 ， 


率 法 线 流 算 子 对 每 个 点 进行 


这 样 可 以 更 加 


有 效 保留 树木 旨 


以 避免 在 分 枝 处 产生 畸形 的 连 
0 在 收缩 点 云 的 同 
， 根 据 现 实 世 界 树木 生物 生 


局 


点 云 时 ， 


收缩 ， 适 当 减 小 每 个 点 的 位 
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数 化 平滑 ， 用 线性 系统 来 控制 几何 平滑 。 在 迭代 平滑 中 ， 通 
过 项 点 加 权 来 实现 全 局 平滑 以 及 保持 网 格 体积 不 变 。 随 着 迁 
代 的 进行 ， 点 云 体积 越 来 越 小 ， 切 平面 的 计算 也 变 得 更 加 的 
困难 ， 并 且 ， 原 始 局 部 网 格 包含 更 多 的 树木 点 云 特征 和 邻 域 
关系 信息 ， 原 始 一 环 邻 域 中 三 角形 较为 工整 ， 可 以 使 整个 收 
缩 过 程 保持 稳定 。 所 以 本 文 只 在 算法 最 开始 的 时 候 计 算 一 次 
邻 域 信 息 ， 得 到 一 环 邻 域 参与 整个 帮 代 过 程 。 
1.1 曲面 曲率 计算 

在 网 格 收缩 的 文献 [10,11] 中 ， 大 多 数 采 用 约束 成 本 最 小 
化 函数 。 下 式 是 最 常见 计算 曲面 8 的 总 曲率 的 函数 : 


CC)=| 尼 +dS (1) 


其 中 : 5 是 网 格 表 面 ，h 和 如 是 曲面 主 曲 率 。 通过 在 顶点 上 拟 
局 部 多 项 式 插值 来 估计 式 (1)。 选择 平滑 交 格 的 一 阶 微分 函 
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oe 


(a) 原 始点 云 (b) 局 部 网 格 化 
(b)local point 


cloud meshing 
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(e) Edge collapse 
图 1 树 点 云 骨 架 提 


Fig.1 Tree point cloud skeleton extraction System 
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部 邻 域 投影 到 二 维 平 面 , 然后 使 用 Delaunay 三 角 剖 分 对 其 进 
以 前 的 网 格 平滑 09 方法 用 所 
通过 求解 拉 普 拉 斯 系统 来 平 
4 格 。Nealen 等 人 07 提 出 一 ee 
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曲率 法 线 流 算 子 


， 真 实 重 现 了 原始 树 


且 对 


提取 的 骨架 进行 骨 


木 点 云 拓扑 结构 和 


缩 。 为 了 将 网 格 平滑 算法 直接 应 
网 格 。 首 先 对 点 云 建 立 KD-tree， 计 算 每 
主 成 分 分 析 算 法 将 局 


如 图 1(b) 所 示 。 


部 网 格 化 和 网 格 3 


F 滑 相关 算法 来 实现 


其 中 :，X 代表 网 格 表面 ，“ 和 ， 为 网 格 上 相互 垂直 的 参数 。 
当 参 数 w=* 时 ， 式 2) 的 函数 值 等 于 总 曲率 。 其 对 应 的 微分 导 
数 为 拉 普 拉 斯 算 子 : 


了 (和 )=X +X， (3) 
其 中 :Xs, 和 X,, 表示 网 格 X 对 4、v 的 二 阶 微分 。 通 过 时 间 t 
对 网 格 X 进行 微分 来 快速 平滑 网 格 高 频 部 分 ， 同 时 保持 网 格 
主要 形状 。 


QX 
F(X) @ 


1.2 隐 式 积分 方案 
前 向 欧 拉 方 案 在 网 格 X 上 对 式 (4) 积 分 : 


Xm =X'+uL(X') (5) 
其 中 :，X2 代表 下 一 次 迭代 的 网 格 ， 和 为 当前 网 格 ， 4 为 时 
间 步 长 。 为 了 满足 算法 的 稳定 性 ， 要 求 x<1， 否 则 ， 网 格 将 
过 度 收缩 ， 会 严重 偏离 网 格 原始 形状 。 


前 向 欧 拉 方 案 思 路 简单 ， 并 且 在 每 个 平滑 阶段 只 需要 线 
性 的 时 间 和 内 存 消耗 。 但 是 当 数 据 集 很 大 时 ， 由 于 时 间 步 长 


的 限制 ， 为 达到 有 效 的 收缩 需要 执行 大 量 的 积分 计算 。 为 了 
节省 时 间 和 解决 数值 稳定 性 问题 ， 本 文采 用 隐 式 积分 方案 ， 
隐 式 积分 是 一 种 稳健 有 效 的 数值 方法 ， 隐 式 欧 拉 方案 有 效 的 
避免 了 时 间 步 长 的 限制 ， 可 以 使 选 代 算法 快速 收敛 。 用 新 生 
成 的 网 格 xm 来 近似 积分 ， 这 样 可 以 更 快 的 达到 平衡 状态 ， 
并 且 具 有 更 好 的 稳定 性 ; 

Xm LX")=X’ (6) 


1.3 构造 曲率 法 线 流 算 子 
在 特定 参数 化 情况 下 ， 
if 算 子 的 近似 值 。 


一 般 采 用 伞 状 算 子 作为 拉 普 拉 
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1 
eR (7) 


: Pj 为 网 格 顶点 ，pi; 为 p; 合 状 邻 域 顶点 ，N() 为 pi 的 


三 网 5 
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在 点 云 记 ， 用 


点 云 


伞 状 邻 域 。 伞 状 介 有 统一 性 优势 ， 如 图 2 所 示 ， 伞 状 算 
村 > 


\ 寞 
是 从 顶点 到 其 邻 域 重心 的 矢量 ,即使 它们 表示 不 同 的 频率 ， 


也 具有 相同 的 拉 普 拉 斯 算 子 值 ,这 表示 它们 具有 相同 的 曲率 ， 
这 意味 着 网 格 应 该 有 具有 相同 长 度 的 演 ， 并 且 顶 点 一 坏 邻 域内 
任意 两 个 相 邻 边 之 间 的 角度 都 相等 ,这 并 不 适用 于 实际 网 格 。 
本 文 需 要 一 个 更 加 精确 的 积分 方程 来 收缩 网 格 项 点 ， 以 
便 仔细 区 分 高 频 和 低频 部 分 ， 保 留 原 始 数据 的 细节 特征 。 为 
此 ， 本 文 设 计 一 个 与 尺度 相关 的 离散 拉 普 拉 斯 算 子 来 更 好 的 
近似 式 (4)。 曲 率 法 线 流 算 子 具有 不 依赖 参数 的 特点 ， 可 以 节 


201904.00013v1 


chinaXiv 


录用 定稿 郝 腾 宇 
省 法 线 计 算 成 本 。 曲 率 法 线 流 通过 沿 着 曲面 法 线 方向 以 平均 
曲率 为 速度 移动 顶点 来 达到 收缩 点 云 的 目的 。 

Op; a 

BE =—kn; (8) 
其 中 :为 点 p; 处 平均 曲率 ，n; 为 pi; 的 法 线 向 量 。 通 常 计算 
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2 ， 树 点 云 骨架 提取 
2.1 We 

ao 等 人 [9 利用 拉 普 拉 斯 算 子 收缩 点 云 提取 一 般 点 云 的 
上 骨架 示 点 云 模 型 的 几何 形状 ， 为 了 使 该 系统 适用 于 


法 向 矢量 的 散 度 来 计算 平均 曲率 ， 如 果 项 点 周围 所 有 面 的 法 
线 都 相同 (该 点 曲率 为 0)， 则 该 项 点 就 不 应 该 再 被 收缩 ， 所 
以 对 曲率 法 线 做 以 下 微分 几何 定义 : 


内 尺 点 二 使 收缩 算法 快速 收敛 ， 本 文 用 尺度 相关 的 曲率 法 
线 流 算 子 社 代替 余 切 算 子 。 如 Au 等 人 0 所 述 ， 通 过 求解 以 
下 线性 系统 来 更 新 顶点 P 的 位 置 


人 


其 中 : 4 点 Pp 周围 最 小 区 域 的 面积 ，V 是 关于 pi; 的 空间 坐标 
(wy,z) 的 导数 。 


Pp: 


(a) 高 频 
(a) High frequency 


(b) 低 频 
(b) Low frequency 
图 2 企 状 算 子 
Fig.2 Umbrella operator 


fa Pp: 
A | \ pg 
Le < 
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(a) 一 环 邻 域 (b) 一 环 邻 域 相 邻 面 
(a) One-ring nelghborhood (b) Adjacent faces in One-ring 


neighborhood 
图 3 曲率 法 线 流 算 子 
Fig.3 Curvature normal flow operator 
为 了 将 式 (9) 应 用 到 网 格 ， 需 要 将 其 离散 化 ， 如 图 3(a) 所 
示 ， 4 为 p; 的 一 环 邻 域 所 有 三 角形 的 面积 。 经 过 从 连续 公式 
版 本 到 离散 版 本 的 推导 09， 得 到 下 式 : 


= Deore teot Blel (10) 

如 图 3(b) 所 示 , % 和 Pb; 是 一 坏 邻 域 中 与 边 % 相对 的 两 个 

角 ，4 是 一 环 邻 域 三 角形 面积 和 。 
本 文 在 原 有 的 拉 普 拉 斯 算 子 的 基础 上 添加 一 环 邻 域 面 积 
盲 息 和 边界 权重 信息 ， 使 得 拉 普 拉 斯 算 子 充分 利用 原始 点 云 


中 : Wi 和 Wi 分别 为 收缩 力 和 吸引 力 的 对 角 加 权 和 矩阵 ，P' 是 
当前 顶点 位 置 ，Pm 是 下 次 更 新 位 置 ， 因 为 上 式 中 只 有 P" 为 
未 知 数 ， 方 程 数 大 于 未 知 数 ， 所 以 可 以 用 最 小 二 乘法 来 解决 
式 (12)。 


长 并 


PerP+ 志 We” 一人 (13) 


其 中 : Wai 表示 Ws 矩阵 第 i 个 对 角 元 素 。 第 一 个 式 子 提供 收 
缩 力 ， 使 用 隐 式 曲率 法 线 流 算 子 沿 顶 点 法 向 方向 除去 局 部 几 
何 细节 特征 ， 第 二 个 式 子 提供 吸引 力 ， 在 迭代 收缩 期 间 保留 
原始 模型 几何 形状 。 


1 


(go) 第 三 次 收缩 


(c) Third contraction 


(a) 第 一 次 收缩 


(a) First contraction 


(b) 第 二 次 收缩 
(b) Second contraction 
图 4 点 云 收缩 
Fig.4 Point cloud contraction 
每 次 迭代 中 ， 基 于 局 部 收缩 比 (初始 单 环 面积 和 当前 单 
面积 比率 》 更 新 每 个 顶点 的 位 置 约 束 Wss 。 如 图 4 所 示 ， 
在 迭代 过 程 中 ， 榈 木 点 云 中 薄 细 的 分 枝 比 粗 密 的 分 枝 率先 达 
到 收敛 条 件 。 如 图 5 所 示 ， 随 着 和 迭代 的 进行 ， 点 云 体积 下 降 
速率 加 快 ， 在 迭代 终止 前 ， 体 积 下 降 速率 逐渐 平缓 。 从 代 
率 法 线 流 算 子 可 以 使 一 环 邻 域内 的 三 角形 保持 良好 的 形状 ， 
从 而 确保 初始 质量 差 的 局 部 单 环 三 角形 正常 收缩 ， 避 免 收 缩 
点 云 发 生 畸 形 。 随 着 迭代 进行 ， 单 环 面 积 越 来 越 小 ， 吸 引力 
权重 越 来 越 大 ， 使 得 系统 矩阵 逐渐 对 角 占 优 。 对 角 占 优 矩 阵 
为 非 奇 异 和 矩阵 ， 所 以 顶点 更 新 过 程 将 越 来 越 稳定 。 在 这 个 过 
程 中 ， 避 免 了 包括 无 穷 大 值 和 除 零 的 数值 误差 。 点 云 收缩 是 


一 


党 


一 个 全 局 约束 的 平滑 过 程 ， 每 次 迭代 都 会 清除 高 频 细节 和 品 


局 部 相关 性 ， 综 合 局 部 领域 信息 突出 单个 点 的 重要 性 ， 反 过 
来 ， 利 用 点 与 点 之 间 的 相关 性 突出 局 域 特征 。 曲 率 法 线 流 算 
子 提供 点 云 收 缩 约束 力 。 通 过 沿 其 近似 曲率 法 线 方向 移动 顶 
点 从 网 格 表 面 去 除 细节 和 噪声 ， 从 而 达到 收缩 点 云 的 目的 。 
求解 离散 曲率 流 算 子 ， 目 的 是 使 点 云 局 部 邻 域 曲率 为 零 ， 即 
LV=0(V 为 项 点 )， 表 示 去 除 顶 点 法 线 分 量 ， 使 项 点 沿 着 其 
法 线 方向 进行 收缩 ,其 中 工 为 曲率 法 线 流 算 子 矩 阵 ，W 为 两 
点 之 间 的 权重 : 


1 
w= (ere +eorp)lel 


了 
Lj= W;= 2 -mw (11) 
(JE 已 


0 


逐渐 变 小 ， 


使 模型 体积 


声 ， 同 时 保留 关键 的 几何 拓扑 特征 ， 
当 体 积 接近 零 时 ， 和 友 代 收敛 。 


图 5 当前 体积 与 原始 体积 的 比例 和 迭代 次 数 的 关系 
Fig.5 Relationship between shrinkage volume ratio and iteration 
2.2 粗 骨 架 提取 
上 一 步 把 点 云 收 缩 为 零 体 积 状 态 ， 接 下 来 将 得 到 的 收缩 
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点 云 C 转 换 为 1D 1 
先 ， 
如 图 


6(b) 所 示 ， 从 


然后 以 5 为 球 心 和 以 7 为 半径 在 C 上 进行 邻 域 查 询 ， 


线 骨架 ， 


Cc 中 随机 取 一 点 G 作 


同时 建立 点 云 -骨架 映射 M 。 首 
] 半 径 为 的 球 对 收缩 点 云 进行 采样 


郝 腾 宇 ， 


》 采样 点 集 为 o 
为 第 一 个 采样 点 3 ， 
球 内 所 


有 点 标记 为 与 采样 点 5 相关 《建立 原始 点 云 与 采样 点 的 映射 


关系 M ),， 在 球 内 


找到 距离 5, 最 远 的 点 5; 为 


依 此 类 推 ， 


关 ， 选 择 距离 C, 近 的 采样 点 标记 。 


直到 所 有 的 收缩 点 被 标记 。 若 C, 与 %& 和 3 同时 相 
采样 完成 后 ， 建 立 原始 点 


云 与 采样 点 的 映射 关系 M ,如 图 6(c) 所 示 , 如 果 C, 与 % 相关 ， 


Cn 与 2 相关 , 则 C, 与 Cn 


共享 一 个 一 环 邻 域 ,连接 5; 和 5;。 若 


存在 孤立 的 采样 点 ， 则 将 它 与 其 最 近 的 采样 点 相连 。 


2 
3 2 二 | 
… 过 ©. 所 村 
Se | | “| 和 
(a) 收 缩 点 云 (b) 球 采样 ”“”(c) 采 样 点 连接 ” (d) 骨 架 网 格 
(a) Shrink (b)Ballsampling (c)Sampling  (d) Skeletal grid 


point cloud 


point connection 


图 6 球 采样 与 采样 点 连 


接 


Fig.6 Ball sampling and sampling point connection 


2.3 ” 折 琶 边 


得 到 粗 骨 架 网 格 之 后 ， 需 要 迭代 折 麦 不 必要 的 边 ， 
骨架 网 格 中 没有 三 角形 。 
点 云 保持 相同 的 拓扑 


为 了 使 提取 的 


架 点 与 输入 点 云 之 | 


结构 ， 同 时 保留 足够 的 骨架 点 来 确保 骨 
司 的 精确 对 应 关系 。 设 


直到 


: 基于 曲率 法 线 流 的 树 点 


下 一 个 采样 点 5;， 


ChinaXiv 合 作 期 刊 


架 提 取 方 法 第 37 卷 第 5 期 


Te(ej)=FCPD)+FCPDD) (16) 

把 @ 折 钱 到 点 p ， 然 后 更 新 p; 的 误差 度量 : 
F(pj)=F(p)+F(p))， 使 之 前 与 p; 相连 的 边 与 p; 相连 。 需 要 注 
意 的 是 , 同时 更 新 原始 点 云 采样 点 映射 M， 将 之 前 与 p; 相关 
的 输入 点 云 更 新 为 与 ?; 相关。 

为 了 避免 骨架 网 格 被 过 度 简化 ， 导 致 生成 的 
过 长 ， 失 去 原始 点 云 与 生成 骨架 的 精确 对 应 关系 。 
成 本 函数 来 惩罚 生成 长 边 的 边 折 车 操 作 。 

LT, (e;)= el lel 


线 骨 架 边 
设置 一 个 


(17) 


如 图 7 所 示 ， 其 中 表示 直接 与 p; 相连 的 点 ， 上 式 表示 
将 p; 合 并 到 pj; 后 ， 与 Pp; 直 接 相连 的 边 的 总 边 长 变化 程度 。 
T(e;)= yal a (ey) + yol, (es) (18) 
7。、%h 分别 为 折合 成 本 了 ole) 和 边 长 成 本 Ts(%y) 的 权重 
占 比 。 原 始 网 格 中 顶点 的 初始 误差 为 0， 当 边 收缩 后 ， 更 新 
顶点 误差 ， 依 次 选取 收缩 后 成 本 最 小 的 边 进行 迭代 收缩 。 总 
成 本 函数 为 形状 成 本 和 抽样 成 本 的 加 权 和 ， 类 似 于 点 云 收缩 
步骤 中 的 收缩 权重 和 位 置 权重 ， 赋 予 形 状 成 本 函数 更 大 的 权 
重 ， 由 于 形状 成 本 是 累加 的 ， 所 以 形状 成 本 逐渐 主导 边 折 芋 
过 程 ， 保 留 骨 架 网 格 的 拓扑 结构 。 
基于 生长 角度 约束 的 骨架 连通 
描 分 辨 率 有 限 以 及 叶子 和 其 他 分 支 的 自我 遮挡 ， 


ly 


HP， 


得 到 的 树木 点 云 会 出 现 断 层 和 不 规则 等 问题 ， 从 而 使 提取 的 


1D 曲线 骨架 与 原始 


个 包含 形状 项 


骨架 不 连通 ， 针 对 该 问题 ， 使 用 双向 探测 方法 连接 主 骨架 和 
其 他 分 支 骨架 。 根 据 树 木 的 生物 特性 ， 不 应 该 存在 角度 变化 
过 大 的 树枝 , 树干 与 分 枝 之 间 连 接 角度 " 不 应 该 超过 45° 。 如 


和 抽样 项 的 成 本 函数 ， 每 次 迭代 折合 最 小 成 本 的 边 。 为 了 节 


约 计算 成 本 ,采用 


7;， 并 且 删 除 包含 边 @ 的 所 有 三 角形 。 


为 了 保留 骨架 


QEM 网 格 简化 方法 作为 形状 成 本 函数 ,通过 计算 每 个 顶点 的 
误差 度量 来 估计 边 折 针 3 


E£ 边 折 奢 ， 折 革 @ 时 将 顶点 pi 合并 到 顶点 
图 7 半边 折 圣 

Fig.7 Half-edge folding 

网 格 的 几何 拓扑 特征 ， 在 网 格 边 上 应 用 


关联 边 2 


间 的 平方 距离 之 和 。 
平方 距离 为 一 CR 
直线 方程 ax+by+cz+4=0 中 的 参数 , 并 日 
为 二 次 基本 误差 矩阵 : 


)p ， 其 中 KE =[apcd 


a ab ac ad 
ab b: bc bd 
dac be ce? cd 
ad bd cd d’ 


p 的 初始 误差 度量 是 pi 与 其 所 有 相 邻 边 


所 以 % 的 误差 


F(pi)= pr 2 (KF K;)p:; 
(beE 


折 车 成 本 为 


起 的 误差 ， 误 差 度量 
点 了 到 其 相关 边 4 所 在 直线 的 


为 顶点 到 其 相 


，(4b,c,d) 为 @ 所 在 


(a +b+e)=1, KIK, 


(14) 


力 的 平方 距离 之 和 : 
(15) 


图 8 所 示 ， 将 含有 最 低 点 的 骨架 定义 为 主 骨架 ， 在 主 骨架 上 
进行 广度 优先 搜索 ， 在 每 个 骨架 点 p; 处 ， 沿 着 从 其 父 节点 到 
该 点 方向 投影 角度 ? 的 圆锥 ， 在 一 定 范围 内 时 ， 当 圆锥 与 子 
骨架 相交 于 点 P ， 将 子 骨架 点 P 拟定 为 与 p; 相连 的 候选 点 。 
在 P 进行 样 条 拟 合 计算 该 点 切线 方向 7， 在 P 处 沿 着 7 投影 

E， 若 圆锥 与 p; 也 相交 ， 则 通过 连接 Pp; 来 连 


团 


一 个 相同 的 


天 


pe 
ss 
(a) 定 位 候选 子 骨 架 (b) 骨 架 连 通 
(a) Positioning candidate sub-skeleton (b) Skeleton connectivity 
图 8 骨架 图 连通 


Fig.8 Skeleton connectivity 

2.5 ”骨架 居中 

在 点 云 收缩 的 迭代 过 程 中 ， 由 于 点 云 密度 不 均匀 ， 导 致 
生成 的 曲线 骨架 偏离 模型 中 心 甚 至 突出 模型 表面 ， 特 别 是 相 
邻 的 部 分 在 密度 和 粗细 程度 上 存在 很 大 差异 的 情况 下 ， 因 为 
密度 和 厚度 较 大 的 区 域 需 要 较 强 的 收缩 力 才 能 收缩 到 零 体 积 
状态 (近似 一 条 直线 ), 这 会 将 其 附近 较 细 的 区 域 的 顶点 拉 向 
较 厚 区 域 的 中 心 ， 从 而 使 薄 细 区 域 的 收缩 点 云 突出 原始 模型 
表 耳 


针对 该 问题， 基于 在 骨架 提取 步骤 中 建立 的 原始 点 云 骨 
架 映 射 M ， 将 每 个 骨架 点 Skeleton 移动 到 其 对 应 的 原始 点 云 
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Mi 的 中 心 ， 如 图 9 所 示 。 从 而 将 每 个 骨架 点 重新 定位 到 其 
局 部 点 云 区 域 的 质心 处 。 
Skeleton, =(Skeleton, + ,12 
eleton, eleton 7 (19) 


\ 


(b) 骨 架 居 中 


(b) Skeleton centered 


(a) 点 云 骨架 映射 
(a) Point cloud-skeleton mapping 
图 9 骨架 居中 
Fig.9 Skeleton centered 

2.6 本 文 算法 流程 

本 文 算法 包括 构建 点 云 一 环 邻 域 . 点 云 收缩 骨架 提取 、 
骨架 居中 、 骨 架 连 通 以 及 几何 重建 。 创 建 一 个 数组 向 量 7 存 
储 点 云 一 环 邻 域 ，Y 为 当前 点 云 体积 ， 初 始 化 收缩 权重 矩阵 
Ws 和 吸引 天 阵风 ， 收 缩 点 云 C ， 曲 率 流 算 子 矩阵 也， 
iter_time 为 迭代 上 限 ，s' 为 p; 一 环 邻 域 面积 ，1 为 当前 迭代 次 
数 , 数 组 Skeleton 存储 骨架 点 ，M 保存 骨架 点 与 原始 点 云 的 映 
射 关 系 ，7 存储 网 格 三 角形 。 

a) 构 建 点 云 一 环 邻 域 ， 计 算 p; 的 KX 近邻 N, ; 对 NN。 进行 
PCA 分 析 和 Delaunay 三 角 训 分 ; 将 pi 一 环 邻 域 存 入 7 ; 

b) 由 式 (11) 计 算 zm; 由 式 (12) 计 算 Cc* ;更 新 收缩 权重 : 
WA = ， 更 新 吸引 权重 : Ws = WS? /Sr ; 


c) 如 果 iter_time <5 或 一 二 >0.01 ， 重 复 步 又 b)， 直至 进 
代 收 敛 ; 

dj) 在 C 上 进行 最 远 点 球 采样 ， 将 采样 点 存 入 Skeleton 中 ， 
根据 邻 域 相关 性 连接 采样 点 ; 
e) 由 式 (18) 更 新 边 权 重 , 折 受 最 小 权重 边 ， 


人 重复 步骤 e)， 直 至 骨架 网 格 中 不 存在 三 角形 ; 
根据 (19) 式 将 骨架 进行 居中 人 处 理 ; 


行 
h) 基 于 本 文 提 出 的 生长 角度 约束 连通 骨架 ; 
3 ”实验 与 分 析 


本 文 的 基于 骨架 提取 的 树木 点 云 重 建 算法 在 CPU i7 2.2 
GHz、 内 存 16 GB 以 及 显卡 Intel Iris Pro 1536 MB 的 Mac 操 
作 系 统 平台 下 实现 。 骨 架 提取 工作 在 Matlab 2017b 平台 上 完 
成 ， 树 木 重建 工作 使 用 Xcode 以 及 OpenGL 库 完 成 。 

相对 于 文献 [18] 的 余 切 算 子 , 图 10 清晰 地 说 明了 曲率 法 
线 流 算 子 的 优势 ， 在 点 云 收 缩 方 面 ， 余 切 算 子 在 细 枝 弯曲 角 
度 较 大 时 ,收缩 点 云 突出 原始 模型 ,并 且 和 忽略 一 些 短小 细 枝 ， 
与 实际 情况 存在 较 大 差异 。 在 提取 的 骨架 方面 ， 细 枝 之 间距 
离 紧密 处 ， 提 取 的 骨架 出 现 环形 结构 现象 ， 另 外 ， 在 分 枝 交 
叉 点 处 ， 存 在 骨架 交叉 问题 。 本 文 算 子 沿 着 点 云 曲 率 法 线 
向 向 点 云 内 部 收缩 ， 即 使 在 树冠 细 枝 处 也 很 好 遵循 该 原则 ， 
不 但 良好 的 保留 树木 的 细 枝 特征 和 几何 形状 ， 而且 避免 了 奇 
异 结构 ， 充 分 表达 树木 的 植物 特性 。 
打破 收缩 力 和 吸引 力 之 间 的 平衡 是 点 云 收缩 的 关键 ， 使 
其 在 总 体 上 提供 一 种 驱动 顶点 向 内 部 收缩 的 力 。 经 实验 证 明 ， 
收缩 尺度 si 直接 影响 点 云 收缩 的 效果 。 如 图 11 所 示 , 若 立 太 
小 ， 点 云 模型 收缩 不 到 零 体积 状态 ， 而 且 需 要 多 次 迭代 和 花 
费 更 多 的 时 间 来 达到 终止 条 件 ; 若 si 太 大 ， 点 云 密 度 大 的 部 
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分 将 吸引 其 附近 密度 小 的 部 分 的 向 其 移动 ， 粗 密 分 枝 和 注 细 
分 枝 交 叉 处 的 收缩 点 云 偏离 模型 中 心 ， 经 过 多 次 实验 验证 ， 
在 提取 树木 点 云 骨架 时 ， 在 [3,5] 得 到 的 收缩 点 云 能 充分 表 
达 模型 的 几何 拓扑 。 


(a) 余 切 算 子 (b) 本 文 算 子 
(a) Cotangent operator (b) Proposed operator 
图 10 算 子 比较 图 
Fig. 10 Operator comparison 


py nMmp WU 


(a) sz =0.7 (b) Y =3 (c) sr =40 

图 11 不 同 大 小 收缩 尺度 的 点 云 收缩 效果 

Fig. 11 Point cloud shrinkage of different contraction Scale 

本 文 算 法 的 重点 是 树木 点 云 骨架 提取 ， 大 多 数 种 类 的 本 
木 的 形状 都 是 细 长 的 ， 本 文 算法 在 4 次 迭代 内 达到 收 勾 。 在 
实验 的 过 程 中 ， 发 现在 对 收缩 点 云 进 行 球 采 样 的 步骤 中 ， 球 
采样 的 半径 的 大 小 直接 影响 提取 的 骨架 的 效果 。 如 图 
12(a)~(c) 所 示 ， 随 着 采样 半径 的 逐渐 增加 ， 提 取 的 骨架 特征 
细节 被 忽略 ， 比 如 ， 树 稍 分 又 处 没有 被 提取 骨架 ， 细 长 分 枝 
变 直 ， 失 去 树木 原 有 的 生物 特性 。 但 是 采样 时 间 和 采样 点 连 
接 时 间 也 相应 减少 ， 另 外 当 树 校 之 间距 离 小 于 采样 半径 时 ， 
产生 环形 结构 。 本 文采 用 自 适应 采样 机 制 ， 根 据点 云 的 密度 
和 分 枝 点 的 数量 自动 设置 采样 半径 ,实现 时 间 和 效果 的 统一 。 
用 户 也 可 以 根据 自己 的 需求 设置 采样 半径 。 实 验证 明 ， 采 样 
半径 r=0.01 时 ， 树 冠 处 细 枝 骨架 清晰 ， 延 伸 方 向 合理 ， 保 留 
有 良好 的 视觉 效果 。 


Wy Ny «tz TY 
pg) WA Ny J 
(a)r=0.01 (br=0.02 (ojr=0.05 


图 12 不 同 球 采 样 半径 的 骨架 对 比 
Fig. 12 The skeletons of different ball sampling radii 
本 文 提 取 的 骨架 充分 保留 了 树木 点 云 的 几何 形状 ， 且 骨 
架 连 通 ， 无 环 状 结构 和 多 余 骨 架 分 枝 。 树 枝 末 段 骨 架 清晰 ， 
树枝 延伸 方向 符合 树木 生长 特性 ， 为 树 的 虚拟 修剪 等 研究 提 
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供 了 数据 支持 。 本 文 算法 不 仅 适应 于 树木 点 云 ， 对 大 多 数 非 
平面 点 云 也 有 相当 显著 的 效果 。 本 文 的 算法 对 于 噪声 点 云 有 
很 好 的 鲁 棒 性 ， 在 点 云 收缩 步 又 中 ， 基 于 隐 式 积分 方案 对 点 
云 表面 进行 平滑 , 用 最 小 二 乘 系统 求解 (14) 式 得 到 收缩 点 云 ， 
时 声 和 细微 的 特征 被 均匀 平滑 ， 使 本 文 算法 对 噪声 不 敏感 ， 
忆 此 ， 不 需要 任何 前 置 滤波 器 对 输入 点 云 进行 预 处 理 。 
图 13 与 14 显示 了 本 文 算法 与 rosa 算法 RI 以 及 Ll 中 轴 
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繁多 及 密度 不 均匀 的 树木 点 云 , Ll 中 轴 算 法 采样 不 均匀 ,于 
失 大 量 主要 枝条 的 特征 信息 ， 骨 架 不 连续 情况 严重 ， 产 生 了 
大 量 的 环形 骨架 ， 枝 条 在 树冠 处 生长 方向 凌乱 无 序 ， 与 现实 
树木 几何 形状 出 入 较 大 , 在 一 般 模型 和 噪声 模型 方面 , L1 提 
取 的 骨架 大 多 数 不 连 通 ， 并 且 在 环形 模型 提取 的 骨架 容易 突 
出 模型 表 函 
相 比 之 下 ， 本 文 算 法 不 需要 任何 预 处 理 和 法 线 信息 


o 


R21 算法 的 比较 结果 。 该 实验 中 ， 本 文 取 f=0.01。rosa 算法 在 
拓扑 结构 复杂 的 树木 点 云 方 面 上 花费 的 时 间 是 本 文 算法 的 几 
十 倍 ， 并 且 提 取 的 骨架 出 现 不 连通 分 枝 ， 树 干 分 又 处 的 骨架 
出 现 畸 形 现象 ， 在 末 枝 处 产生 大 量 奇异 骨架 点 ， 枝 条 骨架 不 
平滑 ， 出 现 枝条 延伸 方向 不 合理 与 相交 连通 现象 ， 在 骨架 末 
尾 处 出 现 分 义 现象 ， 产 生 多 余 骨 架 断 枝 ， 另 外 ，rosa 算法 输 
入 点 云 都 需要 经 过 去 除 噪声 和 异常 值 等 预 处 理 ， 并 且 要 求法 


接 在 原始 点 云 上 进行 操作 ， 节 省 计算 成 本 ,微调 收缩 尺 

采样 半径 就 可 以 避免 产生 拓扑 错误 的 骨架 分 枝 和 环形 骨架 。 
在 不 完整 点 云 处 理 上 ， 点 云 密 度 高 度 不 均匀 ， 提 取 的 骨架 偏 
向 于 点 云 密 度 较 大 的 部 分 ， 偏 离 模 型 中 心 ， 本 文 算法 在 后 期 
工作 中 对 上 骨架 进行 居中 操作 。 在 噪声 点 云 处 理 上 ， 本 文采 取 
隐 式 积分 对 点 云 表 面 进行 平滑 ， 没 有 奇异 骨架 点 产生 ， 本 文 
算法 对 噪声 点 云 有 良好 的 鲁 棒 性 ， 提 取 的 骨架 结构 完整 ， 连 


线 估 计 , 花费 额外 时 间 成 本 。L1 中 轴 算 法 在 原始 点 云 的 下 采 
样 点 上 进行 ， 采 样 点 的 好 坏 直 接 影响 骨架 的 效果 。 对 于 村 条 


(a) 原 始点 云 (bjrosa 算法 
(a) Original point cloud (b) ROSA algorithm 


图 13 不 同 算法 树木 点 云 骨架 提取 结果 对 比 图 


续 、 无 交叉 现象 、 无 环形 结构 、 生 长 方向 合理 、 冠 层 处 骨架 
清晰 ， 真 实 再 现 输 入 模型 的 几何 形态 和 生物 特征 。 


(OL1l 中 轴 算 法 
(c) Ll-axis algorithm 


(d) 本 文 算法 
(d) Procposed algorithm 


Fig. 13 Comparison of different algorithm tree point cloud skeleton extraction results 
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(d) 本 文 算法 
(d) Proposed algorithm 


(OL1l 中 轴 算 法 
(c) Ll-axis algorithm 


图 14 不 同 算法 在 不 同类 型 点 云 骨架 提取 结果 对 比 


Fig. 14 Comparison of different algorithms for extracting results from different types of point cloud skeletons 
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表 1 显示 了 本 文 算法 不 同 数目 的 点 云 在 各 个 阶段 所 用 时 
间 ， 从 表 中 可 以 看 出 ， 在 点 云 数 量 合理 的 情况 下 ， 本 文 算法 
吾 效 率 满足 一 般 用 户 需 求 。 在 构造 一 环 邻 域 阶段 和 点 云 收 
缩 阶段 耗费 大 量 时 间 ， 并 着 点 云 数量 的 增加 ， 点 云 收缩 
阶段 所 用 时 间 增 长 幅度 较 大 ， 所 以 本 文 建议 对 时 间 有 要 求 的 
用 户 ， 可 以 对 输入 点 云 进 行 下 采样 等 预 处 理 。 表 2 显示 了 本 
文 算法 与 L1- 中 轴 算 法 、rosa 算法 的 速率 与 效果 ， 本 文 算法 
不 仅 在 骨架 提取 速度 方面 分 别 是 L1- 中 轴 、rosa 的 3 倍 和 17 
倍 ， 并 且 在 骨架 与 原始 树木 点 云 重 合 度 以 及 树木 枝条 重建 率 
方面 提高 了 25%。 

表 1 各 个 阶段 所 用 时 间 


本 


: 基于 曲率 法 线 流 的 树 点 云 


Table 1 Time at each stage 

点 数 一 环 邻 域 点 云 收缩 边 折 县 总 时 间 

1598 1.483s 0.781s 0.053s 2.317s 
4770 5.057s 3.721s 0.146s 8.924s 
12108 12.676s 24.888s 0.322s 37.886s 
34648 36.718s 71.907s 1.108s 109.733s 
110423 114.349s 1200.146s 26.23s 1340.725s 
1598 1.483s 0.781s 0.053s 2.317s 


表 2 不 同 算法 骨架 提取 的 效率 


Table 2 Efficiency of skeleton from different algorithms 


Ll-Medial rosa 本 文 算法 
提取 速度 0.211k/s 0.034k/s 0.579k/s 
重合 度 0.536 0.731 0.942 
枝条 重建 度 0.473 0.713 0.946 
4 结束 语 
本 文 提出 一 种 基于 es ea te 


算法 ， 本 文 算法 利用 点 云 自 身 信息 ， 直 接 在 点 云 上 操作 ， 


需要 进行 额外 预 处 理 和 法 向 量 计算 ， 对 于 噪声 点 云 和 on 
云 有 良好 的 鲁 棒 性 。 与 之 前 算法 相 比 ， 本 文 算法 在 时 间 方 面 
提高 3 倍 以 上 ， 枝 条 重建 度 提 高 了 25%， 充 分 的 保留 了 树 校 
细节 特征 ， 生 动 形 象 的 还 原 了 真实 树木 形态 ， 不 仅 适 用 于 各 
种 类 型 和 大 小 的 树 ， 也 适用 于 其 他 物种 的 一 些 简单 模型 。 
在 未 来 ， 本 文 算法 需要 加 快 提取 骨架 的 速度 ， 尝 试 提 


带 有 树叶 的 树木 点 云 的 骨架 。 另 外 ， 和 希望 将 点 云 -骨架 映射 
用 到 蒙 皮 动 画 方面 ， 实 现 模拟 树木 在 大 自然 各 种 环境 下 的 
动 和 生长 。 
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